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ABSTRACT 

We investigate the intrinsic shapes of starless cores in the Orion GMC, using the 
prestellar core sample of Nutter & Ward-Thompson (2007), which is based on submil- 
limeter SCUBA data. We employ a maximum-likelihood method to reconstruct the 
intrinsic distribution of ellipsoid axial ratios from observations of the axial ratios of 
projected plane-of-the-sky core ellipses. We find that, independently of the details of 
the assumed functional form of the distribution, there is a strong preference for oblate 
cores of finite thickness. Cores with varying finite degrees of triaxiality are a better fit 
than purely axisymmetric cores although cores close to axisymmetry are not excluded 
by the data. The incidence of prolate starless cores in Orion is found to be very infre- 
quent. We also test the consistency of the observed data with a uniform distribution of 
intrinsic shapes, where oblate and prolate cores of all degrees of triaxiality occur with 
equal probability. Such a distribution is excluded at the 0.1% level. These findings 
have important implications for theories of core formation within molecular clouds. 
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1 INTRODUCTION 

The intrinsic shapes of molecular cloud cores encode im- 
. portant clues about the processes of molecular cloud frag- 
' mentation and core formation. If cores form as the result of 
converging turbulent flows (e.g., see review by MacLow & 
Klessen 2004) then they are expected to have random, triax- 
ial shapes, with a slight preference for prolateness (Gammie 
et al. 2003; Li et al. 2004). If, on the other hand, magnetic 
fields are dynamically important, and ambipolar diffusion 
mediates the core formation process, then the cores are ex- 
pected to be flattened along the magnetic held due to the 
extra magnetic force acting perpendicular to the field lines. 
Core shapes in this case are oblate (Mouschovias 1976) al- 
though not necessarily axisymmetric (Basu & Ciolek 2004; 
Ciolek & Basu 2006). Unfortunately, the distribution of core 
shapes cannot be easily derived from observations, because 
it is not possible to directly deproject the shape of each 
core - statistical techniques have to be used instead. Early 
work on core shapes, assuming axial symmetry, seemed to 
favour prolate cores (Myers et al. 1991; Ryden 1996). How- 
ever, subsequent investigations which relaxed the axisym- 
metry assumption have consistently yielded triaxial, prefer- 
entially oblate core shapes (e.g. Jones et al. 2001; Jones and 
Basu 2002; Goodwin et al. 2002), independently of tracer or 
core sample. 



In light of the recent debate concerning the core for- 
mation timescale and its use as a test for theories of core 
formation in molecular clouds (Hartmann et al. 2001; Tas- 
sis & Mouschovias 2004; Ballesteros-Paredes & Hartmann 
2006; Mouschovias, Tassis & Kunz 2006), the timing is es- 
pecially opportune to revisit the question of the intrinsic 
core shapes, as these may serve as an independent test of 
theories of core formation. The recent progress in studies of 
turbulence-driven core formation and the convergence of the 
findings of different groups concerning the expected distribu- 
tion of core shapes in this scenario (Gammie et al. 2003; Li et 
al. 2004) also offers an unprecedented opportunity for direct 
comparison between theory and observations. We approach 
the problem using a new statistical method: a maximum- 
likelihood analysis, aimed to reconstruct the intrinsic shape 
distribution of molecular cloud cores. Maximum-likelihood 
analyses are particularly powerful because they allow treat- 
ment of different functional forms of the intrinsic shape dis- 
tribution, distinguishing between functional forms, and ex- 
plicitly account for observational uncertainties and different 
assumptions about potential biases in viewing angles. Al- 
though here we present a first treatment under the simplest 
possible assumptions, the expansion of the formalism to in- 
clude more complicated scenarios is straightforward and we 
will pursue it in a future publication. 

We use the recently released core dataset of Nutter & 
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Ward- Thompson (2007) (NWT07), which ofTers the oppor- 
tunity to work with a sample with a series of unique fea- 
tures: a high level of completeness (NWT07 use this sam- 
ple to derive the core mass function in the mass range 0.3 
to ~ IOOM0); cores of varying sizes and masses identified 
in a single survey [in this SCUBA (Submillimeter 

Common User Bolometer Array) survey of the Orion Gi- 
ant Molecular Cloud (CMC)]; confident discrimination be- 
tween prestellar and protostellar cores (protostellar cores 
were identified by the authors using Spitzer data, and we 
have excluded such cores from our analysis). The sample 
used here consists of the 286 prestellar cores from the Orion 
A North and South and Orion B North and South star 
forming regions. Axial ratios for these objects are calculated 
from the quoted semi-major and semi-minor dimensions in 
NWT07. 



2 FORMALISM 

Consider a system of coordinates centered on a triaxial el- 
lipsoid model molecular cloud core with principal axes (a, b, 
c). The triaxial ellipsoid surface in this system obeys 

x'^ 

^ + 12+^ = 1- (1) 
a'' 0^ 

The image of the core appears on the plane of the sky as an 
ellipse, the properties of which can be calculated from Eq. 
IT]) and the orientation of the observer's line of sight (Binney 
1985). If a > fe > c and we define the axial ratios = 6/a 
and ^ = c/a, so 1 > C > ^, the plane-of-the-sky isophotes of 
the core will be coaxial ellipses of axial ratio 
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where 9 and (j) are the line-of-sight orientation angles and 
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The axial ratio q has values limited between and 1. A very 
elongated ellipse has q <^ 1, while a circle has q = 1. 

Let the cores of the Orion CMC have a single intrinsic 
distribution of shapes, porion(C! ?)'^C'^Ci defined as the frac- 
tion of cores with axial ratios between ^ and C + i^C ^ and 
^ -|- d^. Let us assume that porion can be parametrized by 
a set of parameters Cfl From an assumed poiion(C,£,, Ci), 
we can calculate the distribution of axial ratios of projected 
cloud ellipses, p{q), by assuming that the viewing angles are 
random {cf> and cos^ follow uniform distributions). For this 
calculation we use the following Monte-Carlo method: 

(i) We select a set of parameters d for the distribution 

POrion(CC,C'0- 

(ii) We randomly draw a pair of {(, ,^) from the distri- 
bution porion(C) Ci C'Oj a pair of {9, (j>) from a uniform 
probability distribution (equal probability for any value of 



^ For example, in the case of a Gaussian, C\ and C2 would be 
the mean and the standard deviation respectively. 



COS 6 or <j) in the intervals [—1,1] and [0, 27r] respectively). 
We use the values of i^, ^, 9, to calculate q through Eq. ([2} 

(iii) We repeat steps (i)-(ii) a large number of times. We 
bin the obtained values of q, and we normalize the obtained 
distribution p{q) so that it is a probability density function 
(pdf) (/^°°p(g,ft)dg = l). 

Once we have obtained p{q,C'i), we can calculate the 
likelihood function C If we have n observations of core ax- 
ial ratios gi , 52, ...gn then, dropping constant normalization 
factors, the likelihood is 



C{C^,qi,q2,....,qn) =]^p(gj, C",) 



(5) 



The maximum-likelihood parameters of the intrinsic shape 
distribution p((", ^, d) are that set of d which maximize the 
likelihood £ for our given set of observed axial ratios. 

Performing a likelihood analysis requires a selection of 
a functional form for the intrinsic shape distribution p{C,0- 
Lacking any a priori knowledge on the functional form of 
P(Ci ?)i have chosen functional forms based on their sim- 
plicity and appropriateness of their properties (defined in fi- 
nite domains, singly-peaked). We will perform this analysis 
using two distinct families of distributions to test the sensi- 
tivity of our result to the assumed functional form. In the 
following sections, we briefly discuss these functional forms. 



2.1 Beta Distribution of Axial Ratios 

Since both and ^ are defined in finite domains, a natu- 
ral choice is the beta distribution, which is routinely used 
to model events which are constrained within finite value 
intervals. The probability density function of the beta fam- 
ily of distributions is p{x) = x°~^{l — x)''~^ /B{a,b), where 
B{a,b) = r"^(l - tf'^dt is the beta function. The dis- 
tribution is defined in the interval [0,1] and its shape is 
controlled by the non-negative shape parameters a, b. When 
a,b > 1 the distribution is singly-peaked. 

We construct a joint distribution p{(^, ^) assuming that 
follows a beta distribution with domain (0, 1] and ^ follows a 
beta distribution with domain (0, ("]. We start from a double 
beta distribution in variables x and y, which vary between 
and 1, with shape parameters {ax,bx) and {ay,by), respec- 
tively. With the change of variables ^ = x and ^ = Cv [which 
maps the (0, 1] domain of y to the desired (0, C,] domain of 
^], we obtain the joint distribution for i^,^: 



p(C,0 = 



(6) 



B{aa:,ba:)B{ay, by) 

The joint probability of Eq. ((6)1 is such that and ^ are not 
independent. The distribution of Eq. © has a single peak 
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C1C3 —by + 1 

by are such that 
< Co < 1- In Eq. © Ci = (1 - by)/{ay + by - 2), 
C2 = ax — ay — by, and C3 = a^^ ~ Oy + b^ — by — 1. In 
our analysis, we only admit values of the shape parameters 
that result in a singly-peaked distribution. This is a tetra- 
parametric distribution and the likelihood analysis is aimed 

bx. Oy, and by. 
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2.2 Modified Lognormal 



We construct this distribution by seeking an appropriate 
minimal transformation which transforms the domain of 
C, and ^ from (0,1) to (—00,00); we then take the trans- 
formed variables to follow a Gaussian distribution. Such 
an appropriate transformation \8 C, ^ x and £, ^ y where 
X = ln[(^/(l — ()] and y = \n[^/{( — ^)]. Assuming a double 
Gaussian pdf of x and y, the joint pdf of and ^ is 
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1-C) 



(8) 



This is also a tetra-parametric distribution and the like- 
lihood analysis is aimed at determining the values of xq, yo, 
ax, and ay. Due to the transformation used to project the 
(0, 1) interval to (—00, 00), the distribution in ^) space is 
skewed with respect to the original Gaussian shape of the 
distribution in {x, y) space. The distortion increases as a^ 
and ay increase, and eventually the distribution acquires ad- 
ditional local maxima. To avoid such extreme distortions we 
limit the range of ax and ay by imposing a prior requirement 
in our likelihood analysis that ax, ay < 0.7. 



3 RESULTS 

The upper panel of Fig. [1] shows surface plots of the 
maximum-likelihood distribution for the beta-family of dis- 
tributions. Its shape parameters were obtained by maximiz- 
ing the likelihood function, using the "Simulated Annealing" 
algorithm (Corana et al. 1987). On the C ~ C plane we have 
indicated the dividing line between oblate and prolate cores 
(oblate cores have (" > 0.5(1 -I- £,)). On this plane, spheres 
live at = ^ = 1, infinitesimally thin axisymmetric disks 
have — 1, ^ = 0) and infinitesimally thin cigars have 
= 0, ^ — 0). Finite-thickness axisymmetric disks live 
along the C = 1 line. The most likely distribution of intrin- 
sic shapes is peaked at ~ 1 and £, ~ 0.54, indicating that 
the most probable shape is a finite-thickness disk with small 
deviations from axisymmetry. However, larger degrees of tri- 
axiality are also frequently encountered, as the distribution 
decreases only gradually with decreasing (. Even so, only a 
very small fraction of objects in this distribution are prolate. 
The distribution peaks strongly at oblate shapes. 

To test the sensitivity of our likelihood analysis to the 
assumed functional form of the intrinsic shape distribu- 
tion, in the lower panel of Fig. [T] we plot the maximum- 
likelihood modified lognormal distribution. Although the de- 
tailed shape of the two maximum-likelihood distributions is 
different, their overall features are in remarkable agreement. 
The distributions are strongly peaked at about the same 
point (C ~ Ij C ~ 0.5), indicating that the most probable 
intrinsic shape of cores found in the Orion star-forming re- 
gions is predominantly oblate, and the half-thickness of the 
cores is finite (about one-quarter the diameter of the disk). 

Qualitatively this result is to be expected, as most of 
the elongation measurements indicate projected ellipses with 
low elongations, while there are also a few cores with signif- 
icant elongations. If viewed from various random angles, an 
oblate core will most of the time yield an only mildly elon- 
gated projection (face-on or almost face-on observations). 



sphere 




oblate 



sphere 




oblate 



Figure 1. Upper panel: maximum-likelihood modified beta- 
distribution; lower panel: maximum-likelihood modified lognor- 
mal distribution 



while in the few cases when the core is viewed edge-on or 
close to edge-on, a significantly elongated projection will re- 
sult. In contrast, a prolate core viewed at random angles 
will most of the time yield a significantly elongated projec- 
tion, a result not corroborated by the data (see Fig. [3] and 
discussion below). 

Figure[2]gives a sense of the uncertainty associated with 
our determination of the intrinsic shape distribution param- 
eters. It shows contours of the likelihood function for the 
beta-family of distributions. In the upper panel, we plot 
contours of the likelihood as a function of l/ux and 1/bx, 
with Uy and by fixed at their maximum-likelihood values. 
The {ax, bx) parameters primarily control the shape of the 
distribution along the axis. The straight line represents the 
division between the regime where the peak of p((", ,^) cor- 
responds to an oblate object (lower right), and the regime 
where the peak corresponds to a prolate object (upper left), 
assuming that the other two parameters {ay, by) are fixed at 
their maximum-likelihood values. The contours correspond 
to a reduction of the level of the likelihood with respect to its 
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Figure 2. Contours of the likelihood function at a level of 1/3, 
~ 1/22 and 1/100 smaller than the value at the maximum. 
Upper panel: ay, by are fixed at their maximum-likelihood values; 
lower panel: a^, bx are fixed at their maximum-likelihood values. 



maximum value of ^ 1/3, ~ 1/22 and ~ 1/100, which would 
be the levels of the 1,2 and 3a contours for a 2-parameter 
likelihood with log £ falhng as a chi-square distribution. The 
likelihood is very strongly peaked at the oblate part of the 
plane, with the likelihood of a prolate peak being infinitesi- 
mal. 

The lower panel shows the same-level contours for the 
likelihood as a function of 1/ay and 1/by, with ax and bx 
fixed at their maximum-likelihood values. The ay, by param- 
eters primarily control the shape of the distribution along 
the ^ axis. In this case, there is no dividing line since the 
distribution peak would correspond to an oblate object for 
any pair of {ay, by) for the given values of ax,bx: the lat- 
ter fix the peak along the ( axis very close to = 1, and 
objects with such high ( are always oblate independently of 
the value of ^. The likelihood is concentrated close to the 
diagonal, and has a preference for comparable values for ay 
and by, indicating that, along the ^-axis, the intrinsic shape 
distribution is roughly symmetric. The slight ringing at the 
lower-left corner of the plot is a numerical effect due to the 
finite resolution in our calculation of the likelihood, and does 
not affect the values of our maximum-likelihood parameters 
since it occurs sufficiently far from the likelihood maximum. 

Figure |3] shows the distribution of axial ratios for 
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Figure 3. Distribution of axial ratios (minor/major axis) of pro- 
jected starless core ellipses (datapoints). The lines show the pro- 
jected core ellipse distributions derived from a uniform intrin- 
sic shape distribution (solid line), the maximum-likelihood beta- 
distribution (dotted line) and the maximum-likelihood lognormal 
(dashed line). 



the projected core ellipses derived from three intrin- 
sic shape distributions: the maximum-likelihood beta- 
distribution (dotted line), the maximum-likelihood lognor- 
mal distribution (dashed line), and a uniform distribution, 
POrion(CiC) =const. (solid line). The crosses correspond to 
the NWT07 data, binned in intervals of Aq = 0.05. The 
uniform intrinsic shape distribution roughly approximates a 
distribution of core shapes expected from turbulence-driven 
cloud fragmentation and core formation (e.g. Gammie et al. 
2003; Li et al. 2004). That the maximum-likelihood distri- 
butions are a better fit to the data is immediately obvi- 
ous. We have also performed a Kolmogorov-Smirnov test to 
quantify the discrepancy between the data and the uniform 
distribution, and to verify that the maximum-likelihood dis- 
tributions are in fact likely in absolute terms to reproduce 
the observed dataset. We find that the probability that the 
observed data come from the uniform intrinsic shape distri- 
bution is smaller than 0.1%. The same test indicates that 
the dataset is consistent within the la level with originating 
from either of the maximum-likelihood distributions, which 
are peaked at oblate intrinsic shapes. 

A potential source of bias in this analysis is the res- 
olution limitations of the data. It is conceivable that, for 
smaller-sized cores, the actual semi-major and semi-minor 
axes of the projected core ellipse are both comparable to, or 
smaller than, the smallest angle that can be resolved in the 
data. Such a case would lead to an observational bias toward 
circular core ellipses for the smallest objects, in which case 
our statistical analysis would lead to erroneous results. We 
have checked that this is not the case in our dataset. The 
largest number of our smallest elongation cores were not, 
in fact, close to the resolution limit. In addition, as we can 
see in Fig. [3] there is no "pileup" of cores towards g = 1, as 
would be expected if the bias discussed here were significant. 
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4 DISCUSSION 

We have derived the maximum-hkehhood beta-distribution 
of intrinsic shapes for the Orion GMC starless molecular 
cloud cores of the NWT07 set. The maximum- likelihood dis- 
tribution peaks at the shape of an oblate, finite-thickness, 
nearly axisymmetric disk. Oblate cores of varying thickness 
and varying degrees of triax;iality also occur frequently, how- 
ever the distribution falls very quickly to zero in the prolate 
part of the spectrum of shapes. Cores in Orion are most 
likely intrinsically oblate. This main result is robust with 
respect to differences in the assumed functional form of the 
intrinsic distribution. Our result is consistent with the find- 
ings of Jones et al. (2001) and Jones & Basu (2002), who 
used different datasets and a different statistical treatment. 
Hence, it appears that the result that molecular cloud cores 
arc preferentially oblate is also robust to the details of the 
dataset used, the wavelength of the observations, and the 
details of the statistical treatment. A uniform distribution 
of shapes with equal numbers of oblate and prolate cores is 
rejected by the data at a high confidence level (< 0.1%). 

These results are particularly interesting since they can 
have important implications for theories of core formation. 
Turbulence-induced core formation tends to favor uniform 
triaxial shape distributions with, if anything, a slight pref- 
erence for prolate starless cores (e.g. Gammie et al. 2003; Li 
et al. 2004) . Such distributions are not preferred by the ob- 
servational data. In contrast, magnetically-driven fragmen- 
tation naturally leads to oblate objects due to the extra 
magnetic support perpendicular to the field lines. 

An additional effect beyond the core formation process 
which could, in principle, alter the shape of molecular cloud 
cores and introduce a bias favoring oblate shapes is rotation. 
However, the ratios of rotational and gravitational energy 
in prestellax cores such as the ones examined here are ob- 
served to be very small (typically < few percent) and there 
are no observational indications of rotational flattening (e.g. 
Goodman ct al. 1993; Casclli et al. 2002). Even those small 
amounts of angular momentum can lead to appreciable flat- 
tening and rotationally supported discs after sufficient con- 
traction under angular rnorncntum conservation - however, 
such flattening becomes important in much later contraction 
stages, and at small (< 10^ AU) scales (overall, protostel- 
lar core shapes are not significantly different than those of 
prestellar cores, e.g. Goodwin et al. 2002). Similar results on 
the absence of rotationally-induced flattening are also found 
in simulations (Gammie et al 2003; Li et al. 2004). 

The maximum-likelihood approach used in this work is 
a powerful statistical technique, which can be appropriately 
expanded to explicitly account for observational uncertain- 
ties and be combined with magnetic field orientations from 
polarization data as those become available. 

Other than observational uncertainties, our results may 
be biased by three additional factors. First, in order to de- 
rive shapes from the signal-to-noise submillimeter maps, one 
makes the implicit assumption of core isothermality, so that 
emissivity contours truly correspond to the spatial distribu- 
tion of core mass. This source of uncertainty is not of great 
concern, however, because on the one hand isothermality is 
likely to be an excellent approximation for starless cores es- 
pecially at the outer edges of the core defining its shape. On 
the other hand, deviations from isothermality would only 



affect the shape determination if the temperature gradient 
does not follow the density gradient, which is not expected 
for gravitational heating. A second potential complication is 
that these cores all originate in the same GMC. If there is 
a preferential direction in this cloud set by the large-scale 
magnetic field and if the magnetic field is dynamically im- 
portant, then cores may be preferentially oriented with their 
smallest ellipsoid axes aligned with the magnetic field, which 
might weaken our assumption of random viewing angle ori- 
entations. However, in the turbulent core formation scenario, 
this is not a concern since the magnetic field is in this case 
dynamically unimportant and tangled. Therefore, the com- 
parison of the observed core distribution at least with the 
predictions of the turbulence theory of core formation is not 
affected by such a bias, and this effect cannot ameliorate 
the discrepancy between roughly uniform shape distribu- 
tions and the observed data. Additionally, the survey area 
from which our dataset was derived spans a 10° region in 
the sky, rather than a single small cloud where orientation 
biases would be expected to be most significant. Finally, it 
is conceivable that the intrinsic shape distribution is not 
singly-peaked, but involves two distinct peaks in different 
regions of the (C, ^) parameter space. 
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